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ABSTRACT 

For pure fourth order (£ oc R^) quantum cosmology the Wheeler-DeWitt 
equation is solved exactly for the closed homogeneous and isotropic model. 
It is shown that by imposing as boundary condition that \1/ = at the origin 
of the universe the wave functions behave as suggested by Vilenkin. 

PACS numbers: 0450,0460 



1 Introduction 

Non-linear modifications of tlie Eintein-Hilbert action, liave a long history 
[1], and they are of interest, among others, for the following reasons: first, be- 
cause the effective gravitational action predicted by closed bosonic, heterotic 
and supersymmetric strings give the usual Einstein term plus a correction 
quadratic in the curvatures [2]. Second, these theories can be renormilized 
when quantized [3]. Third, pure gravity inflationary models emerge on adding 
an R^ term to the usual gravitational Lagrangian [4]. On the other hand, 
some non-linear Lagrangians can be chosen with the property that the field 
equations for the metric are second order, these are the so-called Lovelock 
actions [5] which can be regarded as formed by the dimensional continuation 
of the Euler characteristics of lower dimensions [6,7,8]. These Gauss-Bonnet 
terms seem to be of importance for the quantization of these theories. 

In this paper we consider a toy quantum model for the theories mentioned. 
The action will be constructed only with the square of the Ricci scalar and the 
corresponding Wheeler-De Witt equation for the closed isotropic cosmologi- 
cal model will be solved by separation of variables. The Vilenkin conditions 
will result after choosing the vanishing of the wave function at the beginning 



of the universe, as first proposed by DeWitt [9]. Previous work considering 
non-linear gravitational Lagrangians has been made by Hawking and Lut- 
trell [10], Boulware [11], Buchbinder and Lyanhovich [12], and Vilenkin [13] 
among others. 



2 Wheeler DeWitt Equation 

The R + 7/?^ theory of gravity, at the classical level is related to Einstein's 
theory of gravity with a scalar field as the source of gravity by a conformal 
transformation [14], a similar relation exists for the scalar-tensor theories of 
gravitation [15]. This relation has been exploited by Kasper [16] to obtain 
approximate solutions to the pure quantum cosmology of fourth order. 

The Theory studied by Kasper for the spatially closed Robertson- Walker 
metric. 



ds^ = -c^dr + a\t) 



+ r [dB + sm "-'-^ ^ 



(1) 



1 — r^ 

where a{t) is the expansion factor of the universe, leads to the following form 
of the Wheeler -De Witt equation in fourth order quantum cosmology: 



12e6"+3'3 + i44e^"+2/3] ^(a, /?) = 0. (2) 
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where a = lna,/3 = Ini? and R is the scalar curvature; therefore this model 
is restricted to i? > 0, since the above equation was obtained using the 
signature -+++, classical anti deSitter space is excluded. Notice that the 
lapse function in that paper is taken equal to 1 in contrast to the paper of 
Hawking and Luttrell where it is proportional to the scale factor. Therefore 
the time coordinates are different in the two cases. 

In order to solve the above equation it is useful to introduce the following 
change of independent variables, 



X 



g2a+/3 _ ^2^^ y = a + (3 = ln(ai?), (3) 



the new variables extend over the half plane a; > (since i? > 0). The 
resulting differential equation is now 



l^i& + aS^ + |: + '2^'-"*^-l*'^'^> = °^ « 



We look for a solution by separation of variables, 



ij{x,y) = F{x)G{y). (5) 

After substitution into the differential equation we obtain 



G' = uG, (6) 



xF" + {u + 1)F' + [-144a; + 12x^]F = 0, (7) 

where u is the separation constant. 

The solution for G is trivial, G = G^e'^y. The solution for F is a little 
more complicated, ii u ^ ±1 it can be obtained as an infinite series. But, if 
z/ = ±1 it is possible to obtain closed form solutions. 

2.1 u = —1 Solution 

In this case Eq.(|^) reduces to 



F" + [-144 + 12x]F = (8) 

This equation can be transformed into a Bessel equation of order 1/3 in the 
variable u = -12^/^ + 12i/3x, F = P{u), 



P" + mP = 0. (9) 



The solution can be written as follows 



F = FiAi{-u) + F2Bi{-u) (10) 

here Ai and Bi are the Airy functions and Fi,F2 are constants that will 
be fixed by the boundary conditions. The complete state function of the 
universe is 



i{j{x, y) = e~y[ciAi{-u) + C2Bi{-u)] (11) 

Taking the limit of this solution when x is large results into the WKB solution 
obtained in reference [16] 



F -. (-^)-V4,±|.3/^ = (-(12i/3x)-i/^e±i(-i2^''-)^'' (12) 

The particular solution of this section is also important when looking at 
the asymptotic solution {x — > oo) in the case of arbitrary value for z/ as will 
be shown in the next section. 



2.2 Boundary condition and WKB approximation 

If we restore the original variables a and R in Eq. (|lT]) we have 



ip[a,K) = — , [i6) 



at a = the values of the Airy functions in the numerator of the wave 
function of the universe are finite in contrast with the denominator that 
vanishes. Therefore the only sensible thing to do is choosing the constants ci 
and C2 so that the whole numerator vanishes at a = 0, and that corresponds 
to the following choice, 



With this choice the numerator in Eq.( |1^) dominates the behaviour at the 
beginning of the universe and the wave function vanishes at a = 0. 

We recall here that the vanishing of the wave function of the universe 
at a = is the boundary condition which Bryce DeWitt [9] argued must 
be imposed on it, because it has the effect of keeping the wavepackets away 
from the singularity. On the other hand if we consider Eq.(D and apply the 
WKB method to it, taking into account the right boundary conditions for 
the case where the potential goes to infinity, that is the case here because 
our coordinate x is defined only for positive values, we obtain. 



S^nh[r,^{l2{l2-x))dx\ 
^{^)WKB = [12(12 -x)]V4 ' (15) 



Now the WKB solution vanishes at a = (x=0). The behaviour of 
the wave function that we have obtained is similar to the hydrogenic wave 
functions when the angular momentum / 7^ in the sense that they vanish 
at the origin of the radial coordinate. 



2.3 u = ISolution 

If in Eq.(|^) we eliminate the first derivative term, to put the equation into 
its normal form by means of the transformation 



F{x)=x-'-'+^^/^gix). (16) 

the differential equation is now 



/+[_12(l2-x) + ^^]^ = 0, (17) 



and we notice that for z/ = 1 the differential equation for g(x) is the same as 
Eq.dl). Therefore the complete wave function of the universe is 



X 



Here we also take as the boundary condition the vanishing of the wave func- 
tion at a = 0, that means that the constants q are given by Eq. ( |14|) . 
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3 Series Solutions 

Before looking at the series solution to Eq. (6) for z^ 7^ ±1, it is useful to 
consider the asymptotic limit x ^ 00 . For that purpose, we look once more 



at the normal form of the differential equation |17 



1 -/y2 
u" + [-12(12 - x) + ^^]u = (19) 



and it is clear that in the asymptotic region [x -^ cxd) this equation reduces 
to the equation solved exactly in the previous section. 

In the case u y^ ±1 we can use Frobenius method assuming a solution for 
Eq.(6) of the form 



F{x) = |a;|" Y. ««^"- (20) 

n=0 



The indicial equation is given by 



a{a + iy) = 0, (21) 

and the solutions are 
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ai = 0,a2 = -i^. (22) 



The recurrence relations are 



a„(n + a — ai){n + a — 02) = 144a„_2 — 12a„_3 (23) 

after substituting the solutions to the indicial equation we have 



a„(n + a){n + a + u) = lUan-2 - 12a„_3. (24) 

From the above equations the series for both solutions is completely deter- 
mined in case u is not an integer; if that is the case the second solution can 
be obtained by the standard procedure used in Frobenius method. 



4 Final remarks 

The simplified model used in this paper has allowed us to exactly solve, for 
the first time, the quantum cosmology of fourth order for a closed Friedmann 
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model. The analogous problem for spatially flat model was solved by Reuter 
[17]. Actions close related to Chern-Simmons are of particular interest at 
the present time; we consider our present work as an attempt in the program 
of studying the quantum cosmology of these more general Lagrangians of 
topological origin. 
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